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ABSTRACT 


New  and  recently  developed  concepts  and  Ideas  useful 
in  obtaining  efficient  computer  algorithms  for  solving  the 
equations  of  motion  of  multi-body  mechanical  systems  with 
flexible  links  are  presented  and  discussed.  These  ideas 
Include  the  use  of  Euler  parameters,  Lagrange's  form  of 
d'Alembert's  principle,  generalized  speeds,  quasi-coordinates, 
relative  coordinates,  structural  analysis  techniques  and  body 
connection  arrays.  The  mechanical  systems  considered  are 
linked  bodies  forming  a  tree  structure,  but  with  no  "closed 
loops"  permitted.  An  explicit  formulation  of  the  equations 
of  motion  is  presented. 


INTRODUCTION 


This  report  discusses  the  development  of  new  methods  for  including  the 
effects  of  flexibility  and  link  and  joint  compliance  in  the  governing  dynami¬ 
cal  equations  of  multi-body  systems.  Specifically,  new,  computationally- 
oriented  techniques  with  potential  for  efficient,  automated,  and  comprehensive 
analyses  of  multi-body  system  dynamics  are  presented  and  discussed. 

The  development  of  equations  of  motion  of  multi-body  mechanical  systems 
has  received  considerable  attention  of  analysts  in  recent  years.  There  are 
several  reasons  for  this:  Foremost,  is  the  fact  that  many  mechanical  systems 
and  devices  can  be  effectively  modelled  by  systems  of  linked  bodies.  But, 
another  reason  is  the  fact  that  it  has  just  recently  been  possible  with  the  aid 
of  high-speed  digital  computers,  to  obtain  efficient  numerical  solutions  of  the 
governing  dynamical  equations.  Hence,  the  emphasis  of  researchers  and  analysts 
working  with  multi-body  systems  has  been  the  formulation  of  equations  of  motion 
which  can  easily  be  developed  into  numerical  algorithms  for  a  computer  code. 

Most  of  the  recent  efforts  in  obtaining  these  dynamic  formulations  and  their 
corresponding  computational  algorithms  has  been  with  systems  of  linked  rigid 
bodies.  Recently,  however,  a  few  researchers  have  attempted  to  include  the 
effects  of  flexibility,  compliance,  and  relative  translation  of  the  links  by 
using  a  variety  of  approaches  such  as  quasi-static  methods,  finite-element  meth¬ 
ods,  modal  analysis,  and  the  strategic  positioning  of  the  flexible  bodies, (for 
example,  to  the  extremities  of  the  system).  Many  of  these  efforts  and  the 
corresponding  methodologies  have  been  stimulated  and  motivated  by  specific 
application  areas  such  as  mechanism  vibration  and  flexible  satellite  oscillations. 
In  this  report,  these  ideas  are  used  and  extended  in  the  outline  of  new  proce¬ 
dures  for  efficiently  modelling  the  dynamics  of  multi-body  systems  with  flexible 
links  and  joints. 
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If  a  mechanical  system  consists  of  connected  bodies  such  that  no  closed 
loops  or  circuits  are  formed,  the  system  is  called  a  "general-chain",  "open- 
chain",  or  "open-tree"  system.  Figure  1.  depicts  such  a  system.  References 
r 1-81  ]*  provide  summary  of  approaches  taken  to  obtain  efficient,  computer- 
oriented  formulation  of  equations  of  motion  for  such  systems  and  related  systems. 
If  the  mechanical  system  model  of  Figure  1.  is  generalized  to  include  transla¬ 
tion  and  compliance  at  the  joints,  it  might  appear  as  shown  in  Figure  2. 
References  2-143  provide  a  summary  of  approaches  taken  to  Include  the  effects 
of  flexibility,  translation,  and  link  and  joint  compliance  of  these  systems. 

In  one  of  these  approaches  it  is  shown  [28,  29,  37,  38,  39,  51]  that  it  is 
possible  to  obtain  expressions  for  the  governing  equations  in  a  form  where  the 
coefficients  are  directly  related  through  computer  algorithms.  This  approach 
uses  Lagrange's  form  of  d'Alembert's  principle,  as  exposited  by  Kane  and  others 
[30,  52,  144,  145,  146]  together  with  body  connection  arrays  [37,  38,  39]  and 
relative  orientation  coordinates  [31,  32,  37,  50]  to  obtain  the  governing  equa¬ 
tions.  Lagrange’s  form  of  d'Alembert's  principle  -  a  virtual  work  type  approach 
combines  the  computational  advantages  of  both  Newton's  laws  and  Lagrange's  equa¬ 
tions.  That  is,  it  has  the  advantage  of  automatic  elimination  of  non-working 
internal  constraint  forces  but  without  the  introduction  of  tedious  differentia¬ 
tion  or  other  similar  calculations. 

Recently,  it  has  been  suggested  by  Huston,  et.  al.  [32,  33,  37]  that  further 
efficiencies  could  be  obtained  through  the  use  of  Euler  parameters  as  described 
by  Whittaker  [ 146]  and  Kane  and  Likins  [ 147] ,  together  with  the  quasi-coordinates 
suggested  by  Kane  and  Wang  [148].  Specifically,  it  is  claimed  [32,  33,  37]  that 
using  Euler  parameters  together  with  relative  angular  velocity  components  as 
generalized  coordinate  derivatives  allows  for  the  avoidance  of  geometrical  singu- 

*  Numbers  .in  brackets  refer  to  References  at  the  end  of  the  report. 
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larities  encountered  with  using  Euler  angles  or  dextral  orientation  angles  to 
define  the  relative  orientation  of  the  bodies. 

The  use  of  Lagrange's  form  of  d'Alembert's  principle,  body  connection  arrays, 
relative  orientation  coordinates,  quasi-coordinates,  and  Euler  parameters  also 

promises  to  provide  an  effective  and  efficient  approach  in  the  modelling,  gover¬ 
ning  equation  formulation,  and  analysis  of  multi-body  systems  with  flexible  links 
and  joints.  The  exposition  of  these  ideas  is  a  primary  objective  of  this  report. 

The  balance  of  the  report  is  divided  into  five  parts  with  the  first  two 
parts  containing  the  geometrical  and  kinematical  development.  The  governing 
dynamical  equations  are  developed  in  the  third  part.  This  is  followed  by  an 
analysis  of  the  flexibility  and  compliance  effects  in  the  fourth  part.  The 
final  part  contains  a  discussion  of  the  developed  procedure  together  with  con¬ 
cluding  remarks. 


PRELIMINARY  GEOMETRICAL  CONSIDERATIONS 


Body  Connection  Array 

Consider  a  mechanical  system  such  as  depicted  in  Figure  1.  To 
develop  an  accounting  routine  for  the  system's  geometry  arbitrarily  select 
one  of  the  bodies  as  a  reference  body  and  call  it  B^.  Next,  number  or 
label  the  other  bodies  of  the  system  in  ascending  progression  away  from 
as  shown  in  Figure  1.  Now,  although  this  numbering  procedure  does  not 
lead  to  a  unique  labeling  of  the  bodies,  it  can  nevertheless  be  used  to 
describe  the  chain  structure  or  topology  through  the  "body  connection 
array"  as  follows:  Let  L(k) ,  k-1, . . . ,N  be  an  array  of  the  adjoining 
lower  numbered  body  of  body  B^.  For  example,  for  the  system  shown  in 
Figure  1.,  L(k)  is: 

L(k)  -  (0,1, 1,3, 1,5, 6, 7, 6)  (1) 

where 

.  (k)  -  (1,2, 3, 4, 5, 6, 7, 8, 9)  (2) 

and  where  0  refers  to  an  inertial  reference  frame  R.  It  is  not  difficult 
to  see  that;  given  L(k) ,  one  could  readily  describe  the  topology  of  the  system. 
That  is.  Figure  1.  could  be  drawn  by  simply  knowing  L(k).  It  is  shown  in 
the  sequel  that  L(k)  is  useful  in  the  development  of  expressions  of 
klnematical  quantities  needed  for  analysis  of  the  system's  dynamics. 
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Transformation  Matrices 


Next,  consider  a  typical  pair  of  adjoining  bodies  such  as  and 
as  shown  in  Figure  3.  The  general  orientation  of  B^  relative  to  B^ 
may  be  defined  in  terms  of  the  relative  orientation  of  the  dextral 
orthogonal  unit  vector  sets  n^  and  (i»l,2,3)  fixed  in  B^  and  Bfc 
as  shown  in  Figure  2.  Specifically  n^  and  n^  are  related  to  each 
other  as 


Bj  i  SJKim5km 


where  SJK  is  a  3  x  3  orthogonal  transformation  matrix  defined  as  [47]: 


“*1,  ‘  2ji  ’  Sta 


(Regarding  notation,  the  J  and  K  in  SJK  and  the  first  subscripts  on  the 
unit  vectors  refer  to  bodies  B^  and  B^,  and  repeated  indices,  such  as  the 
m,  in  Equation  (3)  signify  a  sum  over  the  range  (eg.  1,...,3)  of  that 
index.  Thus,  with  a  computer  SJK^  would  be  the  array  SJK(X^l).) 

From  Equation  (3),  it  is  easily  seen  that  with  three  bodies  B^ ,  B^, 
B.,  the  transformation  matrix  obeys  the  following  chain  and  identity  rules 


SJL  -  SJK  SKL 
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and 


SJJ  -  I  -  SJK  SKJ  »  SJK  SJK_1  (6) 

where  I  Is  the  identity  matrix. 

These  expressions  allow  for  the  transformation  of  components  of 
vectors  referred  to  one  body  of  the  system  into  components  referred 
to  any  other  body  of  the  system  and,  in  particular,  to  the  inertial 
reference  frame,  R.  For  example,  if  a  typical  vector,  V,  is  expressed 
as 


V  -  V 


00 


2ki 


(7) 


then 


V±<°>  -  S0K±j  (8) 

where  0  refers  to  the  inertial  frame,  R. 

Since  these  transformation  matrices  play  a  central  role  throughout 
the  analysis,  it  is  helpful  to  also  have  an  algorithm  for  their  derivative, 
especially  the  derivative  of  SOK.  Using  Equation  (3),  and  noting  that  nQ^ 
are  fixed  in  R,  the  following  is  obtained: 

d(S0Ki:J)/dt  -  n0i  .  Rd  n^/dt  (9) 
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where  the  R  in 


d  n  /dt  indicates  that  the  derivative  is  computed  in 
R.  However,  since  the  n^  are  fixed  in  B^,  their  derivatives  may  be 
written  as  x  n^  where  u>fe  is  the  angular  velocity  of  in  R. 
Equation  (9)  may  then  be  written  as: 


d(S0Kij)/dt  -  _eimn,Jkn  ?0m  *  2kj 


(10) 


or  as 


d(SOK) /dt  -  WOK  SOK 

where  WOK  is  a  matrix  defined  as 


(ID 


WOK 


im 


"imn^kn 


(12) 


and  where  are  the  components  of  referred  to  n^  and  e^^  is  the 
standard  permutation  symbol  [ 150] .  (WOK  is  simply  the  matrix  whose  dual 
vector  [150]  is  u^.)  Equation  (11)  thus  shows  that  the  transformation  matrix 
derivative  may  be  computed  by  a  simple  matrix  multiplication. 


Euler  Parameters 


Finally,  consider  describing  the  relative  orientation  of  B.  and  B. 

J 

by  using  the  so-called  Euler  parameters  as  discussed  by  Whittaker  [  147]  and  Kane 
and  Likins  [148]  It  is  well  known  [147]  that  may  be  brought  into  any 
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general  oriencaeion  relative  to  by  means  of  a  single  rotation  about  an 
appropriate  axis.  If  X^  is  a  unit  vector  along  this  axis  and  if  ^  is  the 
rotation  angle,  the  four  Euler  parameters  describing  the  orientation  of 
relative  to  B^  may  be  defined  as: 

eki  “  xki  sin(V2) 

ek2  "  Xk2  sin<V2) 

(13) 

£k3  "  Xk3  aln(8k/2> 

ek4  "  cos  (6k/2) 


where  the  Xfci  (1*1, 2, 3)  are  the  components  of  Xfc  referred  to  n^,  the  unit 
vectors  fixed  in  B^.  Clearly,  the  (i-1,2,3,4)  are  not  independent  since: 


2  2  2  2 
ekl  +  ek2  +  ek3  +  Ek4 


(14) 


These  parameters  may  be  related  to  angular  velocity  components  by 
using  the  transformation  matrices  as  follows:  It  is  shown  in  [147,  148] 
that  SJK  may  be  expressed  in  terms  of  these  parameters  as: 
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SJK  - 


£kl*'eit2“Ek3+€k4  2  (ekl£k2-ek3ek4)  2  ^kl^"**2^ 


2(ekl£k2+€k3ek4)  "£kl+£k2"ek3+€k4  2(ek2ek3“efclek4) 


2(£kl£k3“€k2£k4)  2(£k2£k3+£kl£k4) 


(15) 


Now,  by  solving  Equations  (11)  and  (12)  for  the  angular  velocity  components, 
one  obtains: 

uikl  -  sok21  sok31  +  sok22  sok32  +  sok23  sok33 


<1^2  «  SOK31  SOK1;l  +  SOK32  SOK^2  +  SOK33  SOK13 


(16) 


oik3  -  SOK11  SOK2l  +  SOK^  2  SOK22  +  SOK^  SOK23 

where  the  dot  designates  time  differentiation.  By  using  Equation  (15), 
these  expressions  may  be  used  to  express  the  components  of  the  angular 
velocity  of  relative  to  B^  in  terms  of  the  Euler  parameters  as: 

\l  “  2(£k4  £kl  “  £k3  £k2  +  £k2  "  £kl  Ek4J 


*\2  "  2(£k3  £kl  +  ek4  £k2  “  £kl  Gk3  “  £k2  £k4) 


(17) 


“k3  "  2(‘£k2  Ekl  +  £kl  Ek2  +  ek4  ek3  “  Ek3  ek4) 
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(Regarding  notation,  in  the  sequel  "hats"  refer  to  relative  angular 
velocity  vectors  or  their  components.  That  is  the  represent  the 
angular  velocity  of  in  R  and  represent  the  angular  velocity  of 
relative  to  B^,  its  adjoining  lower  numbered  body.)  Equation  (17) 
may  now  be  solved  for  the  (i*l,...,4)  in  terms  of  the  leading 
to  the  expressions: 


*kl  “  %(ek4  +  ek3  (^2  “  ek2 

Kl  '  ^<-6k3  “kl  +  ek4  “k2  +  Ekl  \3> 

(18) 

*ek3  "  «ek2  <\l  “  Ekl  ^2  +  £k4 
*k4  "  ^(-£kl  “kl  “  ek2  “k2  “  ek3  “k3} 


This  solution  is  quickly  obtained  by  observing  that  if  Equation  (14)  is 
differentiated  and  placed  with  Equation  (17),  the  resulting  set  of  equations 
could  be  written  in  the  matrix  form: 


v  - 

r 

— 

• 

“kl 

ek4 

■*k3 

Ek2 

-Ekl 

Ekl 

“k2 

£k3 

Ek4 

-Ekl 

_ek2 

Ek2 

-  2 

• 

“k3 

”ek2 

Ekl 

ek4 

"Ek3 

Ek3 

>4. 

_Ekl 

ek2 

ek3 

Ek4 

Ek4 

(19) 
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KINEMATICS 


Coordinates 


A  multibody  system  of  N  bodies,  with  translation  permitted  between 
the  bodies  will.  In  general,  have  6N  degrees  of  freedom.  Let  these  be 
described  by  6N  generalized  coordinates  x^  (lm 1,...,6N)  and  let  the  first 
3N  of  these  be  divided  Into  N  triplets  describing  the  relative 
orientation  of  the  successive  bodies  of  the  system.  Let  the 
remaining  3N  x^  also  be  divided  into  N  triplets  representing  the  relative 
displacement  of  the  successive  bodies  of  the  system.  As  before,  let  B^ 
be  a  typical  body  of  the  system  and  let  be  its  adjacent  lower  numbered 
body,  as  in  Figure  3.  The  angular  velocity  of  relative  to  (that  is, 
the  relative  rate  of  change  of  orientation)  may  then  be  written  as: 


“k  "  “Id  Sji 


+  “k2  2j2  +  “k3  2j3 


(20) 


where  n,.  (j*l,...,N.  i-1,2,3)  are  mutually  perpendicular  dextral  unit 
-J1  J 

vectors  fixed  in  B^.  Next,  let  these  bodies  be  displaced  relative  to 
each  other  with  the  displacement  measured  by  the  vector  £.  as  shown  in 
Figure  4.,  where  0^  and  0^  are  arbitrarily  selected  reference  points  of 
Bj  and  B^.  ,  which  is  fixed  in  B^ ,  is  the  connection  point  or  "origin" 

of  B^.  Then  may  be  written  in  the  form: 


15 


Sk  “  ?kl  2jl  +  St2  2j2  +  *k3  ~j3 

Following  Kane  and  Wang  [149]  introduce  6N  parameters  y^  (A-l 
defined  as: 

yi  “  *1  1  6N 

where  the  first  3N  of  these  are  called  "generalized  speeds",  are 

y3k-2  *  “kl 
y3k-l  *  “k2 
y3k  ““k3 

and  the  remaining  3N  are: 

73(W+k)-2  “  Sd 

♦ 

y30*+k)-l  ”  5k2 

♦ 
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In  general.  Equations  (23)  are  non-integrable.  That  is,  they  cannot 
be  Integrated  to  obtain  generalized  orientation  coordinates  x^_2  » 

,  Xjk  .  Thus,  explicit  parameters  80(1  *3^  do  not 

in  general  exist — hence,  the  name  "quasi-coordinates”.  However,  since 
parameters  are  needed  to  relate  the  relative  orientation  of  the  bodies 
to  the  respective  relative  angular  velocities,  let  the  Euler  parameters 
Introduced  in  the  foregoing  section  be  used  for  this  purpose.  Hence,  if  the 
orientation  of  a  typical  body  relative  to  is  described  by  the  four 
parameters  (i-1, . . . ,4) ,  the  geometry  and  kinematics  of  the  entire  system 

may  be  expressed  in  terms  of  the  4N  Euler  parameters  (k-l,...,N;  1*1,..., 4), 
the  3N  relative  angular  velocity  components  (k-l,...,N;  i-1, 2, 3),  and  the 
3N  displacement  components  5^  (k-l,...,N;  i-1,2,3). 


Lar  Velocity 


The  angular  velocity  of  a  typical  body  in  the  inertial  frame  R  is 
readily  obtained  by  the  addition  formula  as  [143] 


A  A 

+  . , .  +  Wj,  (25) 

where  the  relative  angular  velocities  on  the  right  side  of  this  expression 
are  each  with  respect  to  the  respective  adjacent  lower  numbered  bodies  and 
where  the  sum  is  taken  over  the  bodies  of  the  chain  from  B^  outward  through 
the  branch  containing  B^.  The  L(k)  array  introduced  in  the  foregoing  section 
can  be  useful  in  computing  this  sum:  Consider  for  example,  the  system  shown 
in  Figure  1.  The  angular  velocity  of  3^  is: 
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(26) 


«9  "  §1  +  “5  +  §6  +  ?9 

The  subscript  indices  (ie.  9, 6, 5,1)  may  be  obtained  £rom  L(k)  as 
follows:  Consider  L(k)  as  a  function  mapping  the  (k)  array  (See  Equation  (2)) 
into  the  L(k)  array.  Then,  using  the  notation  that  L®(k)  »  (k), 

L^k)  -  L(L°(k) ,  L2 (k)  -  L(L1(k),  Lj  (k)  -  L(Lj_1(k)),  it  is  seen 

(see  Equation  (1))  that: 

L°(9)  -  9,  LX(9)  -  6,  L2 (9)  -  5,  L3 (9)  -  1  (27) 


Therefore,  u0  may  be  written  as: 

*  p 

Up  -  l  u  ,  q  -  Lp(9) 

pto  ~q 


(28) 


Hence,  in  general,  the  angular  velocity  of  may  be  written  as: 

tt*.  -  I  q  -  LP(k)  (29) 

p-0  ~q 

where  r  is  the  index  such  that  Lr(k)  ■  1  and  it  is  obtained  by  comparing 
LP(k)  to  1.  The  index  r  represents  the  number  of  bodies  from  to  B^  in 
that  branch  of  the  chain  system  For  example,  for  the  system  of  Figure  1. , 
if  k-9,  r»3.  Equation  (29)  is  thus  an  algorithm  for  determining  once 
and  L(k)  are  known. 
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By  examining  Equations  (20,  (23),  and  (25)  it  Is  seen  that  may  be 
written  in  the  form 


Cl) 

-  k 


"kirn  y4  -om 


(30) 


where  there  is  a  sum  over  the  repeated  indices  and  where  (k*l, . . . ,N; 

£**1,...,3N;  m-1,2,3)  form  a  block  array  of  coefficients  needed  to  express 

u>.  in  terms  of  n  .  In  view  of  Equations  (3),  (16),  (20),  and  (23), 

—it  -om 

it  is  seen  that  the  elements  of  the  w  ^  array  may  be  obtained  from  the  SOK 
transformation  matrices.  Moreover,  it  can  be  shown  that  the  matching  between 
the  elements  of  the  «  ^  and  SOK  arrays  is  solely  dependent  upon  the  body 
connection  array  L(k). 


To  see  this,  consider  for  example  the  angular  velocity  of  of  the 
system  of  Figure  1:  From  Equation  (25),  is 

co  -  u  +  0)  +  <u  (31) 

~4  -1  ~3  -4 


where  from  Equations  (3),  (20),  and  (23) 

Sl  *  *1501  +  y2  502  +  y3  5o3  ■  h  5.j  So. 

§3  ■  y7  2ll  +  H  Sl2  +  y9  Sl3  *  y6+j  S01.h  So. 

Sl  ■  y10  ?31  *  711  232  +  y12  533  ■  y9+j  s03»)  5o. 


may  be  written  as: 

(32) 

(33) 

(34) 
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Hence,  the  are: 

5m*  4  “I*2*3 

0  i  -  4,5,6 

“*£.  “  S01mA-9  1  "  7’8'9  ““  1*2.3  05) 

S03mi-9  1  *  10»11»12 
0  i  >  12 


where  6  are  the  identity  matrix  components  [150]. 

Next,  consider  that  the  results  such  as  Equation  (35)  may  be  obtained  for 

the  entire  system  of  Figure  1.  or  Figure  2.  from  a  table  such  as  Table  1. , 

where  the  "m"  entries  of  the  array  are  the  column  of  the  transformation 

matrices.  Finally,  note  that  the  non-zero  entries  in  a  typical  row, 

say  the  kttl  row  of  Table  1.  are  obtained  as  follows:  Let  P  »  L(k). 

Then  SOP  is  placed  in  the  k^  column  of  triplets  of  x^.  Next,  let  Q-L(P). 

til  • 

The  SOQ  is  placed  in  the  P  column  to  triplets  of  x^,  etc.  That  is,  SOM 
is  placed  in  column  L^  3(k)  where  M  »  L3(k),  j-l,...,r+l  with  r  determined 
from  Lr(k)  -  1. 

Finally,  it  is  interesting  to  note  that  the  elements  of  the  array 

(and  hence,  the  transformation  matrix  columns  of  Table  1.)  are  components 
of  the  "partial  angular  velocity  vectors"  as  originally  defined  by  Kane  [144], 
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Angular  Acceleration 


The  angular  acceleration  of  B^  in  R  may  be  obtained  by  differentiating 
Equation  (30).  Noting  that  the  nQffl  are  constant,  this  leads  to: 


2k  ^UikJlm  \  +  UkAm  Horn 


(36) 


A  table  containing  the  aj^p^  can  be  constructed  directly  form  the  corresponding 
table  for  the  to^p^.  For  example,  for  the  system  of  Figure  1.,  such  a 
table  is  shown  in  Table  2. 


Mass  Center  Velocities 

The  velocity  and  acceleration  of  the  mass  center  of  a  typical  body 
(k«l,...,N)  may  be  obtained  as  follows:  Let  r^  locate  relative  0^ 
as  shown  in  Figure  4.  Since  0^  is  located  relative  to  Qfc  by  and  if  is 
located  relative  to  0.  by  the  vector  q.  (See  Figure  4.),  then  by  continuing 

J  -•K 

this  procedure,  G^  may  ultimately  be  located  relative  to  a  fixed  point  0  in 
R,  Che  inertial  reference  frame.  For  example,  for  Body  Bg  of  Figure  2.,  the 
position  vector  Fg  of  Gg  relative  to  0  is: 


~8  "  Ii  +  +  Is  +  36  +  h  +  37  +  &  +  38  +  Is  +  He 


(37) 


In  general,  for  Body  B^,  the  position  vector  of  relative  to  0  is: 
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7,8,9  10,11,12  13,14,15  16,17,18  19,20,21  22,23,24  25,26,27 


24 


(38) 
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where  WK^  is  defined  as: 


^ph  "  ”empi  ^Ai  S0Kph 


v^(3W-iOm  ■  <k»l,...,N;  £.-l,...,3N;  a“l,2,3)  .  (43) 

The  elements  of  the  Vj»  array  are  components  of  the  "partial  velocity  vectors 
as  originally  defined  by  Kane  [144]. 

Mass  Center  Accelerations 


Similarly,  by  differentiation  of  Equations  (40) ,  the  acceleration 
of  G,  in  R  is 


i  •  (iu.  h  +  vki»  V5c 


where  the  non-zero  v^^  are,  by  Equations  (41)  to  (43), 


'kirn  “  ^mhi  rVch  +  l  [  ^mh2,  (Csh  +  qgh>  +  WiW5shl 
q-0 


? 3N,  tn-1,  2,  3) 


Wh*r*  19  5 


«.hl  -  -*mpi  (Vi  S0Sh  +  Vl  s6Kp'n> 


EQUATIONS  OF  MOTION 


Consider  again  a  general  chain  system  such  as  shown  In  Figure  2.,  and 
Imagine  the  system  to  be  subjected  to  an  externally  applied  force  field. 

Let  the  force  field  on  a  typical  body  be  replaced  by  an  equivalent 
force  field  consisting  of  a  single  force  F.  ,  passing  through  G,  together 
with  a  couple  with  torque  M^.  Then  Lagrange's  form  of  d'Alembert's  principle 
leads  to  governing  dynamical  equations  of  motion  of  the  form  [38] : 


F2  +  F2*  “  0  2-l,...,6N  (48) 

F^  (2.-1,..., 6N)  is  called  the  generalized  active  force  and  Is  given 
by: 


F2  "  vk2m  Fkm  +  ^klm  I*km 


(49) 


where  there  is  a  sum  from  1  to  N  on  k  and  from  1  to  3  on  m,  and  where  F 


mk 


F2* 


and  are  the  components  of  F^  and  M^  with  respect  to  n^ 

(2-1,..., 6N)  is  called  the  generalized  inertia  force  and  is  given  by: 


V  -  vk2»  Fk£  +  “k2m 


(50) 
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I 


where  the  Indices  follow  the  same  rules  as  In  Equation  (48) ,  and  where 
F^*  and  are  n^  components  of  inertia  forces,  Ffc*,  and  inertia  torques, 
given  by  [145]. 


Ik* 


A  tk 


(no  sum) 


(51) 


and 


2k*  ■  -Jk  ‘  “k  '  2k  x  ’■hi  *  Sic’  (n0  *") 


(52) 


where  m^  is  the  mass  of  and  1^  is  the  inertia  dyadic  of 
relative  to  G^  (k-1, . . . ,H) .  (F*,  with  line  of  action  passing  through 

together  with  M*  are  equivalent  to  the  Inertia  forces  on  [145]. 
Through  use  of  the  shifter  transformation  matrices,  1^  may  be  written 
in  the  form: 


Ik  *  *kmn  ~om  -on 


(53) 


By  substituting  Equations  (36)  and  (44)  into  Equations  (51)  and  (52) 
and  ultimately  into  Equation  (47),  the  equations  of  motion  may  be  written 
in  the  form: 
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*Jtp 


(54) 


where  there  Is  a  sum  from  1  to  6N  on  p  and  where  a.  and  £ 

*P 

given  by: 


are 


aip  'V  Vkpm  Vk2m  +  *kmn  ^kpm  C‘Jkin 


(55) 


and 


+  ®k  VkAm  Vkqm  +  *kam  \jLrn  *\qn 
+  enmh  *kmr  ^qn  ^ksr  ^kin  ^q 

where  there  is  a  sum  £rom  1  to  N  on  k,  from  1  to  6N  on  q  and  s,  and  from 

X  to  3  on  the  other  repeated  Indices. 

Recall  that  the  first  3N  y^  are  relative  angular  velocity  components. 

These  may  be  related  to  the  Euler  parameters  by  N  sets  of  first  order 

equations  of  the  fora  of  Equations  (18). 

Equations  (54),  (20),  and  the  4N  equations  of  the  form  of  Equations 

(18)  form  a  set  of  13N  simultaneous  first-order  differential  equations  for 

the  6N  y^,  the  3N  and  the  4N  Euler  parameters  (h*l,...,N; 

i«*l,...,4).  Since  the  coefficients  a^  and  f^  in  Equations  (54)  are 

algebraic  functions  of  the  physical  parameters  and  the  four  block  arrays 
•  # 

tukAm*  ^klm*  vk£m  comPuter  algorithms  can  be  written  for  the 

numerical  development  of  these  governing  equations.  Moreover,  once  these 
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arrays  are  developed,  the  system  of  equations  consisting  of  Equations  (54), 
(20),  and  4N  equations  of  the  form  of  Equations  (18),  may  also  be  solved 
numerically  by  using  one  of  the  standard  numerical  Integration  routines 


and  a  linear  equation  solver. 

The  development  of  these  computer  algorithms  and  the  numerical  development 
of  Equations  (54)  might  proceed  a3  follows:  First,  let  the  body  connection 
array  L(k)  (See  Equation  (1))  together  with  ♦■he  geometrical  and  physical 
parameters  rfe,  1^,  and  m^  (See  Equations  (38),  (51),  and  (52).)  and 

the  applied  forces  and  moments  and  (See  Equation  (48).)  be  read  into 

/ 

the  computer.  (Let  r 3^  and,  if  desired,  Ffc  and  be  expressed  in 
terms  of  .)  Next,  from  assumed  initial  values  of  form  the 
transformation  matrix  arrays  SOK  using  Equations  (15)  and  (5).  Use  these 
arrays  to  express  rfc,  1^  and  possibly  F^  and  in  terms  of  n^.  Next, 
using  L(k)  and  SOK  write  an  algorithm,  with  Tables  1.  and  2.  as  a  guide,  to 
form  and  .  For  example,  to  obtain  the  non-zero  ialrgm,  observe  that 

if  L(k)  -  p,  then  “  SOP^  (m-1,2,3;  i-3p+l,  3p+2,  3p+3).  Then,  if 
L(p)  ■  q»  L2(k)  -  q  and  Uy9m  m  SOQ^  (m-1,2,3;  i-3q+l,  3q+2,  3q+3). 

This  assignment  procedure  is  continued  until  unity  is  reached  or  r  times 
where  r  is  given  by  LT(k)  -  1  (See  the  remark  following  Equation  (29).). 
v^g  and  v^9m  may  then  be  obtained  using  Equations  (40)  to  (47).  Finally, 
numerical  values  of  the  coefficients  a^  and  f ^  of  the  governing  differential 
equations  (54)  may  then  be  obtained  from  Equations  (55)  and  (56).  These 
equations  may  then  be  integrated  numerically  to  obtain  incremental  values  to  the 
initial  values  of  the  parameters  yp,  e^,  and  x^  (p-1,. . . ,3N+3;  k-l,...,N;  i»l,2,3, 
and  q-1,2,3),  at  the  end  of  a  time  interval,  say  t^.  New  values  of  the 


transformation  matrix  arrays  SOK  may  then  be  obtained  and  the  entire  process 
repeated  until  a  history  of  the  configuration  and  motion  of  the  system  is 
determined . 

The  application  of  these  expressions  and  ideas  in  an  analysis  of  the 
flexibility  and  compliance  effects  is  developed  in  the  following  part  of  the 
report. 


EFFECTS  OF  COMPLIANCE  AND  FLEXIBILITY 


Let  the  tern  "compliance"  refer  to  the  yielding  or  deformation  of  the  system 
due  to  the  externally  applied  forces  and  due  to  the  inertia  forces.  If  the 
assumption  is  made  that  the  compliance  of  a  link  or  joint  is  "small"  compared  with 
the  general  dimensions  of  the  system,  then  the  effects  of  the  compliance  can  be 
determined  directly  from  the  integration  of  Equations  (54). 

To  see  this,  consider  a  typical  integration  step  as  described  at  the  end  of 
the  foregoing  section.  If  the  v^  *k£m>  u*^,  efci,  y£ ,  and  y?  are  known, 

all  the  kinematics  is  known.  That  is,  by  using  Equations  (30)  and  (36)  the  angular 
velocity  and  angular  acceleration  of  each  body  is  determined.  Similarly,  by  using 
Equations  (40)  and  (44),  the  velocities  and  accelerations  of  the  mass  centers  of 
each  of  the  bodies  is  determined.  Then,  by  using  Equations  (51)  and  (52)  the 
equivalent  inertia  force  system  on  each  body  Is  determined.  Hence,  since  the 
externally  applied  force  field  on  each  body  is  also  known,  the  entire  force  system 
on  each  body  Is  determined.  Therefore,  by  taking  successive  free  body  diagrams  of 
the  bodies  of  the  system,  starting  with  the  Nth  body  and  working  backward  through 
the  chain,  the  force  system  transmitted  across  each  connection  joint  may  be  deter¬ 
mined.  Finally,  by  knowing  the  complete  force  system  acting  on  each  body,  includ¬ 
ing  the  forces  transmitted  across  the  connection  joints,  the  physical  force-defor¬ 
mation  relations  may  be  used  to  determine  the  compliance.  Then  by  addition  and 
superposition  the  compliance  of  the  entire  system  is  determined. 

To  illustrate  this  procedure  in  more  detail,  assume,  for  example,  that  the 
bodies  of  the  system  are  long  slender  members  which  can  be  modelled  as  beams  with 
uniform  cross  section.  Hence,  the  system  of  Figures  1.  and  2.  might  appear  as 
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<  ' 


shown  in  Figure  5.,  and  a  typical  member  of  this  system  might  be  depicted  as  in 
Figure  6.,  where  a  rectangular  shape  is  assumed  and  where  an  axes  system  is  intro¬ 
duced.  As  before  0^  is  the  connection  point  of  the  adjacent  lower  numbered  body 
Bj.  and  is  the  connection  point  with  the  adjacent  higher  numbered  body  B^. 

Let  the  forces  exerted  on  Bj  by  B^  and  B^  at  the  connection  joints  Oj  and 

A  A 

be  represented  by  single  forces  t±/j  aIU*  ik/j  pass*n8  through  0^  and 
together  with  couples  with  torques  and  .  Similarly,  let  the  externally 

applied  force  system  on  B^  together  with  the  Inertia  force  system  of  B^  be  repre¬ 
sented  by  equivalent  force  systems  at  the  ends  Oj  and  of  Bj .  Hence,  let  the 
resultant  force  system  exerted  on  B^at  be  represented  by  the  single  force  -k/j 
passing  through  Q^,  together  with  a  couple  with  torque  and  let  f^j  and  n^/j 

be  expressed  in  the  forms: 


and 


f  *  /.  B  f  , i  n , -  f  . .  n._  4*  f  , ,  n 

~k/j  ~jl  -jl  j2  -j2  j3  ~j3 


^k/j  =  “jl  ~jl  +  mJ2  2j2  +  mj3  2j3 


(57) 


(Note,  that  from  equilibrium  considerations,  the  force  system  exerted  on  B^  at  0^ 

is  equivalent  to  a  single  force  -f.  .  passing  through  0.  together  with  a  couple 

-K/j  j 

with  torque  -  m^ ^ .) 

Let  the  displacement  of  relative  to  0^  due  to  the  beam  compliance,  or 

deformation,  be  represented  by  u..  Let  u,  be  written  in  the  form: 

-j  ~j 


Sj  *  UJ1  +  U12  2j2  +  UJ3  Sj3 


(59) 


Similarly,  let  the  rotation  of  the  beam  cross  section  at  relative  to  the  cross 

section  at  0,,  due  to  the  beam  compliance,  be  represented  by  Let  <j>.  be 

3  *0  ~J 
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written  in  the  form: 


ij  *  *jl  2jl  +  *j2  2j2  +  ^j3  2j3 


(60) 


By  following  the  precedures  of  matrix  structural  analysis  [15]  u^  and 


(61) 

(62) 

(63) 

(64) 

(65) 

(66) 


terms 

of  f^  and  m^  (i 

m 

1,  2,  3)  as: 

UJi 

-  U^/AjEj)  fj]L 

“j2 

-  (*j/3 EjJj2)  fj2 

+ 

v-y  “j  v  nj3 

UJ3 

-  (Z*/3 Ej 

- 

•  V  ”3  2 

© 

-  °ji 

(S'j/2EjJj2)  fj2 

+ 

(I3/E3JJ2>  "j3 

*j3 

-  (^/2EjJj3)  fj3 

- 

”j2 

where  is  the  axial  length  of  is  the  cross-sectional  area,  J 


ji 


(1  -  1,2,  3) 


are  the  centroldal  second  moments  of  area  of  the  cross  section  relative  to  the  X 
axes,  Ej  is  the  elastic  modulus,  and  is  the  shear  modulus. 


ji 


u.  and  <P.  thus  represent  the  compliance  as  yielding  of  member  B,  due  to  the 
~J  -J  J 

holding  and  system  motion.  In  an  automated  analysis  as  outlined  in  the  foregoing 
section,  u.  and  4>.  would  be  calculated  and  then  used  to  adjust  the  geometrical 
parameters  at  each  Integration  step.  Specifically,  3^  (See  Figure  4.)  and  SOK  are 
adjusted  as: 

3k  *  3k  +  Sj 


and 


SOK  +S0K  CJ 


(67) 

(68) 


where  CK  is 


C J  - 

"1 

_0j3 

♦32 

*j3 

1 

-♦3i 

5  2 

6n 

1 

(69) 
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(The  development  of  Equation  (69)  follows  from  the  successive  multiplications  of 
matrices  of  the  form  of  Equation  (15)  for  0  being  and  ^j3  (small  angles) 

about  the  X^,  ^j2»  and  ^j3  axes  respectively.)  The  integration  would  then  pro¬ 
ceed  with  the  adjusted  values  of  q,  and  SOK. 

K 

Vibration  and  Impact  Response 

The  above  compliance  analysis  is  a  quasi  static  analysis  and  as  such  it  does 
not  directly  account  for  oscillatory  or  vibration  phenomena  due  to  the  flexibility 
of  the  system  and  the  externally  applied  (for  example,  impact)  and  inertia  forces. 
If,  as  before,  it  is  assumed  that  the  vibrations  have  relatively  small  amplitude, 
then  a  modelling  and  description  of  the  vibration  phenomena  may  be  obtained  through 
torsion  and  translation  springs  introduced  at  the  connection  points,  or  joints,  of 
the  system. 

To  illustrate  this,  consider  again  a  system  consisting  of  long  slender  members 
which  can  be  modelled  as  beams  as  in  Figure  5.  Consider  two  typical  adjoining  mem¬ 
bers  of  such  a  system  as  shown  in  Figure  7.  The  contribution  to  the  systems 
oscillation  due  to  the  flexibility  of  can  be  modelled  by  1)  three  torsion  springs 
connecting  the  surfaces  of  and  0^  with  spring  constants  Gj  5  and  Ej 

and  governing  the  relative  rotation  of  B.  and  B,  about  axes  parallel  to  n. . ,  nJO 

j  k  —  j  J-  ~  j  *- 

and  n._  respectively;  and  by  2)  three  translation  springs  connecting  Q,  and  0.  with 

J  iC  K 

spring  constants  AjE^/£j,  3EjJj3/2,j  and  3EjJ^2/£j  and  governing  the  relative 

translation  of  B.  and  B,  along  axes  parallel  to  n.,,  n.»,  n.„.  (these  constants  are 
j  k  ~jl  -j^  ~j3 

determined  from  elementary  structural  analysis  as  in  Equations  (61)  to  (66).) 


Discussion 

Both  of  the  above  analyses  Involve  the  effects  of  forces  and  moments  trans¬ 
mitted  across  connection  joints.  The  dynamics  analysis  of  the  preceding  part  is 
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particularly  well  suited  for  accomodating  the  Introduction  of  these  forces  and 
moments  and  for  obtaining  their  contributions  to  the  generalized  forces.  That 
lst  although  the  compliance  procedure  above  suggests  the  use  of  successive  free 
body  diagrams  to  obtain  the  force  and  moment  components,  and  although  this  pro¬ 
cedure  could  be  automated,  these  components  as  well  as  the  spring  and  moment  com¬ 
ponents  of  the  above  oscillation  analysis,  can  be  readily  obtained  and  directly 
incorporated  into  the  governing  equations  by  using  the  partial  velocity  and  partial 
angular  velocity  vectors  of  the  preceding  dynamics  analysis. 


To  see  this,  consider  again  two  typical  adjoining  bodies  such  as  and 
of  Figures  4.  and  7.  As  in  Figure  4.,  let  measure  the  displacement  of  0^  rela¬ 
tive  to  Q^.  Then  by  differentiation  of  [1451  the  velocity  of  0^  may  be  ex¬ 
pressed  as:  V°k  -  XlE^  +  n^  (70) 

As  before,  let  the  orientation  and  rotation  of  B^  relative  to  B^  be  defined  in 
terms  of  Euler  parameters  and  relative  angular  velocity  components.  Then,  from 
Equation  (25)  the  angular  velocity  of  may  be  expressed  as: 

2k +“k’Sj  +  “ktSji  <71> 


Let  the  force  system  which  B^  exerts  on  be  equivalent  to  a  single  force 
~k/j  Pa88*n8  through  together  with  a  couple  with  torque  m^^ ,  as  in  Equations 
(57)  and  (58).  Then  by  the  law  of  action-reaction,  the  force  system  which  Bj 
exerts  on  B^  is  equivalent  to  a  single  force  -  Passing  through  0^  together 

with  a  couple  with  torque  -  m^^ 


Let  y^  (£-1,  ...,  6N)  be  the  generalized  coordinate  derivatives  as  defined 
by  Equations  (23)  and  (24).  Let  the  contribution  to  the  generalized  active  force 
Fg  by  these  forces  transmitted  across  the  connecting  joint  be  F^.  Then  F^  is  given 


by  the  expression  [145]: 
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P£  -  (aV^/Sy^)  •  (fk/j)  +  OUj/3^  )  *  (Sv:/j>  (72) 

+02^/3ya)  ‘(-fk/j)  +  <3V3y£)#  (_sk/j) 

Consider  che  following  cases: 

Case  1:  y^  is  not  equal  to  either  £k^  or  In  this  case  the 

partial  velocities  and  partial  angular  velocities  of  Qk» 

0  t  B , ,  and  B,  may  be  expressed  by  using  Equations  (70) 

K  J  K 


and  (71)  as: 

3YQk/3yi  -  3V%&y£ 

(73) 

and 

(74) 

a  yk/a  y£  *  a  ojj  y^ 

A 

Then,  by  Equation  (72),  F£  becomes: 

/V 

F„  _  „ 

(75) 

Case  2:  y^  is  equal  to  one  of  the  (i“l»  2,  3) .  In  this  case, 
the  partial  velocities  and  partial  angular  velocities  of 
Qk>  0k,  Bj,  and  Bk  become: 


-  3fk/9Cki  -  o  (76) 

a -  Sji  (77) 

aojj/ay^  -3^/^Cki  -0  (78) 

and 

3“k  /3yi  “3^Cki  "°  (79) 


Hence,  by  Equation  (72)  and  (57),  becomes: 


Case  3:  is  equal  to  one  of  the  (1*1,  2,  3).  In  this  case,  the 

partial  velocities  and  partial  angular  velocities  of  Q^,  0^, 


Bj,  and  B^  become: 


a/Vsy*  -  -  o 

33^/3  y A  -  3/k/aS^  .  o 


(83, 

(82) 


3~j/37£,  -  *  0 


(83) 


aV3yll  "  *V3\i  "  5ji 


(84) 


Hence,  by  Equation  (72)  and  (58),  becomes: 


Ft  -  5,i  -  -”,i 


(85) 


The  above  three  cases  include  the  contribution  to  the  generalized  active 
forces  for  each  of  the  y^  (2.  -  1,...,  6N) .  Moreover,  each  of  the  non-zero  con¬ 
tributions  (from  Equations  (80)  and  (85)  occurs  individually;  that  is,  each  con¬ 
tribution  occurs  separately  in  one  of  the  governing  equations.  Hence,  if  in  a 
particular  configuration  or  motion  of  the  system,  y^  is  specified  (for  example, 
y^  is  zero)  then  the  %  th  governing  equation  becomes  an  uncoupled  linear  espres- 
sion  for  the  unknown  restraining  force  or  moment,  thus  determining  the  compliance. 
Conversely,  if  y^  is  an  unknown  variable  (representing  a  degree  of  freedom)  the 
contribution  to  due  to  the  flexibility  as  modelled  by  the  translation  and 
torsion  springs  Is  determined  directly  by  Equations  (61)  to  (66) ,  (80) ,  and  (85) . 
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CONCLUDING  REMARKS 


The  results  of  numerically  solving  the  governing  differential  equations  (54) 
where  the  coefficients  are  given  by  Equations  (55)  and  (56)  are  reported  and 
discussed  in  References  [1,  32,  34,  35,  36,  77,  152-155]  for  a  number  of  physical 
systems  and  configurations  (e.g.  human-body  models,  head-neck  models,  and  flexible 

i 

cables) . 


The  application  of  Equations  (54)  with  these  systems,  however,  is  based  on 
the  use  of  relative  orientation  angles  between  the  respective  bodies  of  the  system 
as  the  generalized  coordinates  (x^)  as  opposed  to  the  use  of  Euler  parameters, 
quasi-coordinates,  and  generalized  speeds  as  outlined  herein.  A  problem  which 
arises  In  the  numerical  solution  of  Equations  (54)  where  orientation  angles  are 
used  Is  that  there  always  exists  values  of  the  angles  and  hence,  configurations  of 
the  system,  for  which  the  determinant  of  is  zero.  A  numerical  solution  will, 

of  course,  fail  to  converge  at  these  singular  configurations  of  the  system,  and 
convergence  is  very  slow  for  configurations  in  the  vicinity  of  a  singularity. 

This  problem  is  avoided  by  using  Euler  parameters  to  relate  the  orientation  geom¬ 
etry  to  the  angular  velocity. 


The  advantages  of  using  Lagrange's  form  of  d'Alembert's  principle  to  obtain 
the  governing  equations  of  motion  for  multi-body  mechanical  systems  has  been 
exposited  in  detail  in  References  [28-30].  Basically,  this  principle  has  the 
advantages  of  Lagrange's  equations  or  of  virtual  work  in  that  non-working  internal 
constraint  forces,  between  the  bodies  of  the  system,  are  automatically  eliminated 
from  the  analysis,  and  may  therefore  be  ignored  in  the  formulation  of  the  governing 
equations.  The  principle,  however,  has  the  additional  advantage  of  avoiding  the 
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differentiation  of  scalar  energy  functions.  Indeed,  the  differentiation  required 
to  obtain  velocities  and  accelerations  are  performed  by  vector  cross  products  and 
multiplication  algorithms  —  procedures  which  are  ideally  suited  for  numerical 
computation.  As  with  Lagrange’s  equations,  Lagrange’s  form  of  d'Alembert's  prin¬ 
ciple  requires  the  use  f  generalized  coordinates  to  define  the  system  geometry. 

The  use  of  Euler  parameters  to  avoid  problems  with  singularities,  as  discussed 
above,  leads  naturally  to  the  use  of  generalized  speeds  -  that  is,  relative  angular 
velocity  components  as  the  generalized  coordinate  derivatives.  This,  in  turn, 
leads  to  additional  computational  advantages  as  observed  by  Kane  and  Wang  [149] 
and  Likins  [122].  Specifically,  by  using  generalized  speeds  (relative  angular 
velocity  components)  as  the  principle  parameters  of  the  analysis,  the  coefficient 
matrices  in  the  governing  equations  can  be  obtained  directly  from  the  body  connec¬ 
tion  array  L(k)  (see  Tables  1.  and  2.). 

The  use  of  "relative"  coordinates,  that  is,  angular  velocity  components  of 
the  bodies  with  respect  to  their  adjoining  bodies,  as  opposed  to  "absolute"  coor¬ 
dinates,  (for  example,  angular  velocity  components  in  inertial  space)  also  con¬ 
tributes  to  the  computational  advantage.  In  applications  with  specific  geometrical 
configurations  [1,  31,  32,  34-36,  50,  77,  152-155],  it  is  seen  that  the  geometry  is 
more  easily  described  in  terms  of  relative  coordinates. 

The  generalization  to  allow  translation  between  the  bodies  of  the  system  makes 
the  analysis  applicable  to  a  much  broader  class  of  problems  than  was  possible  with 
those  previous  analyses  which  are  restricted  to  linked  multi-body  systems.  For 
example,  with  the  head-neck  systems  of  References  [152,  154,  155]  the  use  of  trans¬ 
lation  variables  between  the  vertebrae  is  necessary  to  obtain  satisfactory  models 
of  such  systems.  But,  and  of  perhaps  greater  significance,  the  generalization  to 
Include  translation  between  the  bodies  of  the  system  is  necessary  for  an  efficient 
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analysis  of  the  flexibility  and  compliance  effects  as  discussed  earlier. 

In  this  regard,  the  compliance  can  be  modelled  with  a  quasi-static  approach 
whereas  the  oscillations  and  impact  response  require  a  dynamic  analysis  with  the 
Introduction  of  additional  degrees  of  freedom.  In  both  of  these  cases  the  analy¬ 
sis  outlined  herein  (using  Lagrange's  form  of  d'Alembert's  principle  together  with 
the  use  of  generalized  speeds)  accomodates  the  effects  of  flexibility  in  an  extreme¬ 
ly  efficient  manner.  That  is,  the  forces  and  moments  transmitted  across  the  system 
joints  are  directly  determined  and  incorporated  into  the  governing  equations. 
Moreover,  the  modelling  may  be  made  as  detailed  as  necessary  by  introducing  non- 
linearities  through  the  elastic  springs  and  dampers  and  by  increasing  the  number 
of  joints  and  bodies  of  the  system. 

Finally,  the  entire  analysis  outlined  herein  is  developed  with  the  Intent  of 
obtaining  efficiencies  in  a  computer  oriented  development  of  the  governing  dynami¬ 
cal  equations.  As  such,  its  most  productive  application  will  be  with  large  multi¬ 
body  systems  such  as  finite-segment  models  of  the  human  body,  chains,  cables, 
robots,  manipulators,  and  teleoperators. 
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Hew  and  recently  developed  concepts  and  ideas  useful  in  obtaining  effi¬ 
cient  computer  algorithms  for  solving  the  eouations  of  motion  of  multi-body 
mechanical  systems  with  flexible  links  are  presented  and  discussed.  These 
ideas  include  the  use  of  Euler  parameters,  Lagrange’s  form  of  d’Alembert's 
principle,  generalized  sneeds,  ouasi-coordinates,  relative  coordinates,  struc¬ 
tural  analysis  technlaues  and  body  connection  arrays.  The  mechanical  systems 
considered  are  linked  bodies  forming  a  tree  structure,  but  with  no  "Closed 
loooa’^'cermitted.  An  explicit  formulation  of  the  eouations  of  motion  is 
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